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Abstract: In this chapter, we examine several issues for ordering or partially or- 
dering elements of hyper-power sets involved in the DSmT. We will show the benefit 
of some of these issues to obtain a nice and interesting structure of matrix represen- 


tation of belief functions. 


3.1 Introduction to matrix calculus for belief functions 


A’ rightly emphasized recently by Smets in [9], the mathematic of belief functions is often cum- 
bersome because of the many summations symbols and all its subscripts involved in equations. 
This renders equations very difficult to read and understand at first sight and might discourage potential 
readers for their complexity. Actually, this is just an appearance because most of the operations encoun- 
tered in DST with belief functions and basic belief assignments m(.) are just simple linear operations 
and can be easily represented using matrix notation and be handled by elementary matrix calculus. We 
just focus here our presentation on the matrix representation of the relationship between a basic belief 
assignment m/(.) and its associated belief function Bel(.). A nice and more complete presentation of 
matrix calculus for belief functions can be found in [6] [7] [9]. One important aspect for the simplification 
of matrix representation and calculus in DST, concerns the choice of the order of the elements of the 


This chapter is based on a paper [4] presented during the International Conference on Information Fusion, Fusion 2003, 


Cairns, Australia, in July 2003 and is reproduced here with permission of the International Society of Information Fusion. 
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power set 2°. The order of elements of 2° can be chosen arbitrarily actually, and it can be easily seen 
by denoting m the bba vector of size 2” x 1 and Bel its corresponding belief vector of same size, that 
the set of equations Bel(A) = } 2 gc4 m(B) holding for all A C © is strictly equivalent to the following 
general matrix equation 


Bel=BM-m © m=BM'. Bel (3.1) 


where the internal structure of BM depends on the choice of the order for enumerating the elements of 
2°. But it turns out that the simplest ordering based on the enumeration of integers from 0 to 2” — 1 
expressed as n-binary strings with the lower bit on the right (LBR) (where n = |@|) to characterize all 
the elements of power set, is the most efficient solution and best encoding method for matrix calculus 
and for developing efficient algorithms in MatLa it or similar programming languages [9]. By choosing 
the basic increasing binary enumeration (called bibe system), one obtains a very nice recursive algorithm 
on the dimension n of © for computing the matrix BM. The computation of BM for |O| = n is just 
obtained from the iterations up to i +1 = n of the recursive relation [9] starting with BMo = [1] and 


where 0,41 denotes the zero-matrix of size (i + 1) x (i + 1), 





BMi+1 = SDL Sones (3.2) 
BM; BM; 

BM is a binary unimodular matrix (det(BM) = +1). BM is moreover triangular inferior and symmet- 
rical with respect to its antidiagonal. 
Example for © = {01, 02, 63} 
The bibe system gives us the following order for elements of 2° = {ag,..., a7}: 

ao =000=0 ay = 001 = 0: az = 010 = 02 a3 = 011 = 01 U b2 

œ, = 100 = 3 a5 =101=6,U03 &sş = 110 = 2U &z= 111 = ı U 2 U 03 = O 
Each element a; of 2° is a 3-bits string. With this bibe system, one has m = [m(ag),...,m(a7)]’ and 


Bel = [Bel(ao),..., Bel(a7)|’. The expressions of the matrix BM3 and its inverse BM3_' are given by 


100000 0 o0] 
11000000 
10100000 
11110000 
BM; = 
10001000 
11001100 
10101010 


1Matlab is a trademark of The MathWorks, Inc. 
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1 0 0 0 0 0 0 0 

-1 1 0 0 0 0 0 0 

-1 0 1 0 0 0 0 0 

a |1 -1 -1 1 0 0 0 0 
BM;™' = 

-1 0 0 0 1 0 0 0 

1 -1 0 0 -1 1 0 0 

1 0 -1 0 -1 0 1 0 


at 1 1 -1 1 -1 -i 1 
3.2 Ordering elements of hyper-power set for matrix calculus 


As within the DST framework, the order of the elements of D© can be arbitrarily chosen. We denote the 


=| 


Dedekind number or order n as d(n) £ |D®| for n = |O|. We denote also m the gbba vector of size d(n) x 1 
and Bel its corresponding belief vector of the same size. The set of equations Bel(A) = $ gepe pca M(B) 


holding for all A € D® is then strictly equivalent to the following general matrix equation 
Bel=BM-m © m=BM"'!.-Bel (3.3) 


Note the similarity between these relations with the previous ones (3.1). The only difference resides 
in the size of vectors Bel and m and the size of matrix BM and their components. We explore in the 
following sections the possible choices for ordering (or partially ordering) the elements of hyper-power set 
D®, to obtain an interesting matrix structure of BM matrix. Only three issues are examined and briefly 
presented in the sequel. The first method is based on the direct enumeration of elements of DÌ according 
to their recursive generation via the algorithm for generating all isotone Boolean functions presented in 
the previous chapter and in [3]. The second (partial) ordering method is based on the notion of DSm 
cardinality which will be introduced in section[B.2.2] The last and most interesting solution proposed for 
partial ordering over D® is obtained by introducing the notion of intrinsic informational strength s(.) 


associated to each element of hyper-power set. 


3.2.1 Order based on the enumeration of isotone Boolean functions 


We have presented in chapter Ja recursive algorithm based on isotone Boolean functions for generating 
D® with didactic examples. Here is briefly the principle of the method. Let’s consider © = {6;,...,4n} 
satisfying the DSm model and the DSm order u,, of Smarandache’s codification of parts of Venn diagram 
© with n partially overlapping elements 6;, i = 1,...,n. All the elements a; of DÈ can then be obtained 
by the very simple linear equation dn = D, - un where dy = [ao = 9,01,...,Qa(n)-1]' is the vector of 
elements of DÌ, u, is the proper codification vector and D,, a particular binary matrix. The final result 


dn is obtained from the previous matrix product after identifying (+,-) with (U, N) operators, 0-2 with Ø 
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and 1- x with x. Dn is actually a binary matrix corresponding to isotone (i.e. non-decreasing) Boolean 


functions obtained by applying recursively the steps (starting with D§ = [0 1]’) 


e Ds is built from Df, —ı by adjoining to each row r; of Df_, any row r; of D}, such that r;Ur; = rj. 


Then D,, is obtained by removing the first column and the last line of Df. 


We denote r*8°(a;) the position of a; into the column vector d, obtained from the previous enumer- 
ation/generation system. Such a system provides a total order over D® defined Vai, aj € D® asa; < Qj 
(a; precedes a;) if and only if r“*°(a;) < r*8° (aj). Based on this order, the BM matrix involved in (8:3) 
presents unfortunately no particular interesting structure. We have thus to look for better solutions for 


ordering the elements of hyper-power sets. 


3.2.2 Ordering based on the DSm cardinality 


A second possibility for ordering the elements of D® is to (partially) order them by their increasing DSm 


cardinality. 


Definition of the DSm cardinality 


The DSm cardinality of any element A € D®, denoted Cm(A), corresponds to the number of parts of A in 
the Venn diagram of the problem (model M) taking into account the set of integrity constraints (if any), 
i.e. all the possible intersections due to the nature of the elements 6;. This intrinsic cardinality depends 
on the model M. M is the model that contains A which depends on the dimension of Venn diagram, 
(i.e. the number of sets n = |O| under consideration), and on the number of non-empty intersections in 
this diagram. Cm(A) must not be confused with the classical cardinality |A| of a given set A (i.e. the 


number of its distinct elements) - that’s why a new notation is necessary here. 


Some properties of the DSm cardinality 


First note that one has always 1 < Cm(A) < 2” — 1. In the (general) case of the free-model M/ (i.e. the 
DSm model) where all conjunctions are non-empty, one has for intersections: 

Cult) =... = Cms (Ân) = 2”! 

Cms(0i N 8j) = 2"? for n > 2 

Cms lOi N 8j N Ok) = 2"-3 for n > 3 

It can be proven by induction that for 1 < m < n, one has C ms (0i N bia N...96;,,) = 2”7™. For the 
cases n = 1, 2, 3,4, this formula can be checked on the corresponding Venn diagrams. Let’s consider this 
formula true for n sets, and prove it for n + 1 sets (when all intersections/conjunctions are considered 
non-empty). From the Venn diagram of n sets, we can get a Venn diagram with n + 1 sets if one draws 


a closed curve that cuts each of the 2” — 1 parts of the previous diagram (and, as a consequence, divides 
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each part into two disjoint subparts). Therefore, the number of parts of each intersection is doubling 
when passing from a diagram of dimension n to a diagram of dimension n + 1. Q.e.d. 
In the case of the free-model Mf, one has for unions: 

Crys (0i U0;) = 3(2"-?) for n > 2 

Cms (0i U0; U 0k) = 7(2"~3) for n > 3 

It can be proven also by induction that for 1 < m < n, one has Cmr (0i U Oin U... U Oin) = 
(2™ —1)(2"-™). The proof is similar to the previous one, and keeping in mind that passing from a Venn 
diagram of dimension n to a dimension n + 1, each part that forms the union 0; 0; N Ok will be split 


into two disjoint parts, hence the number of parts is doubling. 


For other elements A in DÌ, formed by unions and intersections, the closed form for C qs (A) seems 
more complicated to obtain. But from the generation algorithm of D°, DSm cardinal of a set A from 
D® is exactly equal to the sum of its coefficients in the u, basis, i.e. the sum of its row elements in 
the D, matrix, which is actually very easy to compute by programming. The DSm cardinality plays in 
important role in the definition of the Generalized Pignistic Transformation (GPT) for the construction of 
subjective/pignistic probabilities of elements of DÈ for decision-making at the pignistic level as explained 
in chapter[ZJand in [5]. If one imposes a constraint that a set B from D® is empty, then one suppresses 
the columns corresponding to the parts which compose B in the D, matrix and the row of B and the 
rows of all elements of D? which are subsets of B, getting a new matrix D’,, which represents a new 
model M’. In the u, basis, one similarly suppresses the parts that form B, and now this basis has the 


dimension 2” — 1 — Cy4(B). 


Example of DSm cardinals on M/ 


Consider the 3D case © = {01,02,03} with the free-model Mf corresponding to the following Venn 
diagram (where < i > denotes the part which belongs to 6; only, < ij > denotes the part which belongs 


to 6; and 6; only, etc; this is Smarandache’s codification (see the previous chapter). 


0 02 


Figure 3.1: Venn Diagram for Mf 
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The corresponding partial ordering for elements of DÈ is then summarized in the following table: 


Ae D® Cms (A) 


A 
ao = 0 


en 


ay Ê 61 N b2 N b3 
ay £ 01 N b2 

az 2 61 N b3 

a4 £ b2 N b3 

as £ (01 U 02) N 03 
ag £ (01 U 03) N 02 


a7 £ (02 U 03 NA, 


I> 


) 
) 
) 

ag = (01 N 82) U (81 N 03) U (02 N 03) 

ag = 6; 

10 £ 02 

ai Ê 03 

a2 £ (01 N 02) U 03 

aiz £ (01 N 03) U 62 

aya £ (02 N 03) U Oy 

Q15 £6, Ub 

aig = 6, U b3 


Q17 2 02 U 83 


a åO nn a a on F A A A www N NDN 


Q18 £ 01 U6, U 03 
Table 3.1: Cms (A) for free DSm model M/ 


Note that this partial ordering doesn’t properly catch the intrinsic informational structure/strength 
of elements since by example (81 M62) U (01 N83) U (02 N03) and 6; have the same DSm cardinal although 
they don’t look similar because the part < 1 > in 4; belongs only to 6; but none of the parts of 
(81 N 02) U (01 N 03) U (82 N 43) belongs to only one part of some 6;. A better ordering function is then 
necessary to catch the intrinsic informational structure of elements of D°. This is the purpose of the 


next section. 


Example of DSm cardinals on an hybrid DSm model M 


Consider now the same 3D case with the hybrid DSm model M 4 Mf in which we force all possible 


conjunctions to be empty, but 6; N 02 according to the following Venn diagram. 
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Figure 3.2: Venn Diagram for M 


The corresponding partial ordering for elements of D® is then summarized in the following table: 


Ae D® Cmu(A) 
ao =0 0 
a, ê 0 N bə 1 
ag Ê 63 1 
a3 £90 2 
a4 £ b2 2 
as £ 01 U b2 3 
ag = 01 U b3 3 
az Ê b2 U 03 3 
ag = 81 U 2 U 63 4 





Table 3.2: Cm(A) for hybrid DSm model M 


Another example based on Shafer’s model 


Consider now the same 3D case but including all exclusivity constraints on 0;, i = 1,2,3. This corre- 


sponds to the 3D Shafer’s model M? presented in the following Venn diagram. 


0, b2 


m Ø 


Then, one gets the following list of elements (with their DSm cardinal) for the restricted DÌ, which 


coincides naturally with the classical power set 2°: 
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ag = 0 
ay, = 6; 1 
az Ê bə 1 
az £ 03 1 
as = 01 U b2 2 
as £ 6, U b3 2 
ag = 62 U 03 2 
az Ê 01 U 02 U b3 3 


Table 3.3: Cmo(A) for Shafer’s model M? 


The partial ordering of D? based on DSm cardinality does not provide in general an efficient solution 
to get an interesting structure for the BM matrix involved in (8.3), contrary to the structure obtained by 
Smets in the DST framework as in section [B.I] The partial ordering presented in the sequel will however 


allow us to get such a nice structure for the matrix calculus of belief functions. 


3.2.3 Ordering based on the intrinsic informational content 


As already pointed out, the DSm cardinality is insufficient to catch the intrinsic informational content of 
each element d; of D®. A better approach to obtain this, is based on the following new function s(.), which 
describes the intrinsic information strength of any d; € D®. A previous, but cumbersome, definition of 
s(.) had been proposed in our previous works [I] 2] but it was difficult to handle and questionable with 


respect to the formal equivalent (dual) representation of elements belonging to D®. 


Definition of the s(.) function 


We propose here a better choice for s(.), based on a very simple and natural geometrical interpretation 
of the relationships between the parts of the Venn diagram belonging to each d; € D®. All the values of 


the s(.) function (stored into a vector s) over D® are defined by the following equation: 
8 eq 


s = Dp: wn (3.4) 


with s = [s(do) ... s(dp)]/ where p is the cardinal of DÌ for the model M under consideration. p is 
equal to Dedekind’s number d(n) — 1 if the free-model M/ is chosen for © = {01,..., 0n}. Dn is the 
hyper-power set generating matrix. The components w; of vector Wn are obtained from the components 


of the DSm encoding basis vector u,, as follows (see previous chapter for details about D,, and un) : 


wi Ê 1/U(us) (3.5) 
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where l(u;) is the length of Smarandache’s codification u; of the part of the Venn diagram of the model 


M, i.e the number of symbols involved in the codification. 


For example, if u; =< 123 >, then l(u;) = 3 just because only three symbols 1, 2, and 3 enter in the 


codification u;, thus w; = 1/3. 


From this new DSm ordering function s(.) we can partially order all the elements d; € DÌ? by the 


increasing values of s(.). 


Example of ordering on D°=t%}2} with M/ 


In this simple case, the DSm ordering of D® is given by 


a; E€ D® s(ai) 

ao = 9 s(ao) = 0 

a, =6,N02 s(a,) =1/2 

ag = bı s(a2) = 1 + 1/2 

ag = b2 s(a3) =1+4+1/2 

as =1Ub2 s(a4) =14+1+4+1/2 





Based on this ordering, it can be easily verified that the matrix calculus of the beliefs Bel from m by 


equation (8.3), is equivalent to 


Bel(@) 100 0 0 m(0) 
Bel(61 N 62) 1 1 0 0 0 m(0ı N b2) 
Bel(62) 1 1010 m(02) 
Bel(61 U 62) 1 1 1 1 1 m(0ı U 02) 
x ae N m N RA oY 
Bel BM2 m 


where the BM2 matrix has a interesting structure (triangular inferior and unimodular properties, 


det(BM2) = det(BM;') = 1). Conversely, the calculus of the generalized basic belief assignment 


m from beliefs Bel will be obtained by the inversion of the previous linear system of equations 


58 CHAPTER 3. PARTIAL ORDERING ON HYPER-POWER SETS 


m(0) 1 0 0 0 0 Bel(Q) 
m(61 N 02) -1 1 0 0 0| |Bel(6, N82) 
m(0ı U 02) 0 1 -1 -1 1 Bel(@; U 02) 

x ae N ae N z P2 
m MB2=BM};' Bel 


Example of ordering on D®={1:2283} with Mf 


In this more complicated case, the DSm ordering of D® is now given by 



































a; € D°,i=0,...,18 (aj) 

ø 0 

01 N 02 N 83 1/3 

01.62 1/3+1/2 

01 N 03 1/3 + 1/2 

62 N 03 1/3 + 1/2 

(01 U 02) N 03 1/3 + 1/2 + 1/2 

(01 U 03) N 02 1/3 + 1/2 + 1/2 

(02 U 03) N 01 1/3 + 1/2 + 1/2 

(8, NO2) U (6,9 03) U (82 N03) 1/3 +1/2+1/2+1/2 

6; 1/3+1/2+1/2+1 

b2 1/3+1/2+1/2+1 

b3 1/3+1/2+1/2+1 

(01 N 02) U 03 1/3+1/2+1/2+1+1/2 
(01 N03) U b2 1/3 +1/2+1/2+ 1+ 1/2 
(02 N 03) U 01 1/3+1/2+1/2+1+1/2 
61 U b2 1/3+1/2+1/2+1+1/2+1 
61 U 63 1/3+1/2+1/2+1+1/2+1 
b2 U 63 1/3+1/2+1/2+1+1/2+1 
0, U 62 U b3 1/3+1/2+1/2+1+1/2+1+1 








The order for elements generating the same value of s(.) can be chosen arbitrarily and doesn’t change 
the structure of the matrix BM3 given right after. That’s why only a partial order is possible from s(.). 
It can be verified that BM; holds also the same previous interesting matrix structure properties and that 


det(BM3) = det(BM;') = 1. Similar structure can be shown for problems of higher dimensions (n > 3). 
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Although a nice structure for matrix calculus of belief functions has been obtained in this work, and 
conversely to the recursive construction of BM, in DST framework, a recursive algorithm (on dimension 
n) for the construction of BM, from BM,,_; has not yet be found (if such recursive algorithm exists ...) 


and is still an open difficult problem for further research. 


10000000 00 0 0 0 0 0 0 0 0 0 
1 100000000 0 0 0 0 0 0 0 0 0 
1110000000 0 0 0 0 0 0 0 0 0 
1101000000 0 0 0 0 0 0 0 0 0 
1100100000 0 0 0 0 0 0 0 0 0 
11011100000 0 0 0 0 0 0 0 0 
1110310100000 0 0 0 0 0 0 0 
11110001 00 00 00 0 0 0 0 0 
1 1 1 1 1 1 1 1 1 0 000000000 
BM;=|1 1 1100010100000000 0 
1 1 1 0 10 1000100000000 
1 1 O 1 1 10000010000000 
1 1 1 1 1 1 1 1 1 0 0 1 1000000 
1 1 1 1 1 1 1 1 1 0 1 00100000 
1 1 1 1 1 1 1 1 1 1 000010000 
1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 1000 
1 1 1 1 1 1 1 1 1 1 O 1 10 10 100 
1 1 1 1 1 1 1 1 1 0O 1 1 1 10001 0 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


3.3 Conclusion 


In this chapter, one has analyzed several issues to obtain an interesting matrix representation of the 
belief functions defined in the DSmT. For ordering the elements of hyper-power set D? we propose three 
such orderings: first, using the direct enumeration of isotone Boolean functions, second, based on the 
DSm cardinality, and third, and maybe the most interesting, by introducing the intrinsic informational 
strength function s(.) constructed from the DSm encoding basis. The third order permits to get a nice 
internal structure of the transition matrix BM in order to compute directly and easily by programming 


the belief vector Bel from the basic belief mass vector m and conversely by inversion of matrix BM. 
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